Genotyping, sequencing and analysis of 140,000 adults from Mexico City.

The Mexico City Prospective Study is a prospective cohort of more than 150,000 adults recruited two decades ago from the urban districts of Coyoacán and Iztapalapa in Mexico City1. Here we generated genotype and exome-sequencing data for all individuals and whole-genome sequencing data for 9,950 selected individuals. We describe high levels of relatedness and substantial heterogeneity in ancestry composition across individuals. Most sequenced individuals had admixed Indigenous American, European and African ancestry, with extensive admixture from Indigenous populations in central, southern and southeastern Mexico. Indigenous Mexican segments of the genome had lower levels of coding variation but an excess of homozygous loss-of-function variants compared with segments of African and European origin. We estimated ancestry-specific allele frequencies at 142 million genomic variants, with an effective sample size of 91,856 for Indigenous Mexican ancestry at exome variants, all available through a public browser. Using whole-genome sequencing, we developed an imputation reference panel that outperforms existing panels at common variants in individuals with high proportions of central, southern and southeastern Indigenous Mexican ancestry. Our work illustrates the value of genetic studies in diverse populations and provides foundational imputation and allele frequency resources for future genetic studies in Mexico and in the United States, where the Hispanic/Latino population is predominantly of Mexican descent.

Converging evidence from exome sequencing and common variants implicates target genes for osteoporosis.

In this study, we leveraged the combined evidence of rare coding variants and common alleles to identify therapeutic targets for osteoporosis. We undertook a large-scale multiancestry exome-wide association study for estimated bone mineral density, which showed that the burden of rare coding alleles in 19 genes was associated with estimated bone mineral density (P < 3.6 × 10-7). These genes were highly enriched for a set of known causal genes for osteoporosis (65-fold; P = 2.5 × 10-5). Exome-wide significant genes had 96-fold increased odds of being the top ranked effector gene at a given GWAS locus (P = 1.8 × 10-10). By integrating proteomics Mendelian randomization evidence, we prioritized CD109 (cluster of differentiation 109) as a gene for which heterozygous loss of function is associated with higher bone density. CRISPR-Cas9 editing of CD109 in SaOS-2 osteoblast-like cell lines showed that partial CD109 knockdown led to increased mineralization. This study demonstrates that the convergence of common and rare variants, proteomics and CRISPR can highlight new bone biology to guide therapeutic development.

A deep catalog of protein-coding variation in 985,830 individuals.

Coding variants that have significant impact on function can provide insights into the biology of a gene but are typically rare in the population. Identifying and ascertaining the frequency of such rare variants requires very large sample sizes. Here, we present the largest catalog of human protein-coding variation to date, derived from exome sequencing of 985,830 individuals of diverse ancestry to serve as a rich resource for studying rare coding variants. Individuals of African, Admixed American, East Asian, Middle Eastern, and South Asian ancestry account for 20% of this Exome dataset. Our catalog of variants includes approximately 10.5 million missense (54% novel) and 1.1 million predicted loss-of-function (pLOF) variants (65% novel, 53% observed only once). We identified individuals with rare homozygous pLOF variants in 4,874 genes, and for 1,838 of these this work is the first to document at least one pLOF homozygote. Additional insights from the RGC-ME dataset include 1) improved estimates of selection against heterozygous loss-of-function and identification of 3,459 genes intolerant to loss-of-function, 83 of which were previously assessed as tolerant to loss-of-function and 1,241 that lack disease annotations; 2) identification of regions depleted of missense variation in 457 genes that are tolerant to loss-of-function; 3) functional interpretation for 10,708 variants of unknown or conflicting significance reported in ClinVar as cryptic splice sites using splicing score thresholds based on empirical variant deleteriousness scores derived from RGC-ME; and 4) an observation that approximately 3% of sequenced individuals carry a clinically actionable genetic variant in the ACMG SF 3.1 list of genes. We make this important resource of coding variation available to the public through a variant allele frequency browser. We anticipate that this report and the RGC-ME dataset will serve as a valuable reference for understanding rare coding variation and help advance precision medicine efforts.

Common and rare variant associations with clonal haematopoiesis phenotypes.

Clonal haematopoiesis involves the expansion of certain blood cell lineages and has been associated with ageing and adverse health outcomes1-5. Here we use exome sequence data on 628,388 individuals to identify 40,208 carriers of clonal haematopoiesis of indeterminate potential (CHIP). Using genome-wide and exome-wide association analyses, we identify 24 loci (21 of which are novel) where germline genetic variation influences predisposition to CHIP, including missense variants in the lymphocytic antigen coding gene LY75, which are associated with reduced incidence of CHIP. We also identify novel rare variant associations with clonal haematopoiesis and telomere length. Analysis of 5,041 health traits from the UK Biobank (UKB) found relationships between CHIP and severe COVID-19 outcomes, cardiovascular disease, haematologic traits, malignancy, smoking, obesity, infection and all-cause mortality. Longitudinal and Mendelian randomization analyses revealed that CHIP is associated with solid cancers, including non-melanoma skin cancer and lung cancer, and that CHIP linked to DNMT3A is associated with the subsequent development of myeloid but not lymphoid leukaemias. Additionally, contrary to previous findings from the initial 50,000 UKB exomes6, our results in the full sample do not support a role for IL-6 inhibition in reducing the risk of cardiovascular disease among CHIP carriers. Our findings demonstrate that CHIP represents a complex set of heterogeneous phenotypes with shared and unique germline genetic causes and varied clinical implications.

Multiancestry exome sequencing reveals INHBE mutations associated with favorable fat distribution and protection from diabetes.

Body fat distribution is a major, heritable risk factor for cardiometabolic disease, independent of overall adiposity. Using exome-sequencing in 618,375 individuals (including 160,058 non-Europeans) from the UK, Sweden and Mexico, we identify 16 genes associated with fat distribution at exome-wide significance. We show 6-fold larger effect for fat-distribution associated rare coding variants compared with fine-mapped common alleles, enrichment for genes expressed in adipose tissue and causal genes for partial lipodystrophies, and evidence of sex-dimorphism. We describe an association with favorable fat distribution (p = 1.8 × 10-09), favorable metabolic profile and protection from type 2 diabetes (~28% lower odds; p = 0.004) for heterozygous protein-truncating mutations in INHBE, which encodes a circulating growth factor of the activin family, highly and specifically expressed in hepatocytes. Our results suggest that inhibin ßE is a liver-expressed negative regulator of adipose storage whose blockade may be beneficial in fat distribution-associated metabolic disease.

Germline Mutations in CIDEB and Protection against Liver Disease.

BACKGROUND: Exome sequencing in hundreds of thousands of persons may enable the identification of rare protein-coding genetic variants associated with protection from human diseases like liver cirrhosis, providing a strategy for the discovery of new therapeutic targets. METHODS: We performed a multistage exome sequencing and genetic association analysis to identify genes in which rare protein-coding variants were associated with liver phenotypes. We conducted in vitro experiments to further characterize associations. RESULTS: The multistage analysis involved 542,904 persons with available data on liver aminotransferase levels, 24,944 patients with various types of liver disease, and 490,636 controls without liver disease. We found that rare coding variants in APOB, ABCB4, SLC30A10, and TM6SF2 were associated with increased aminotransferase levels and an increased risk of liver disease. We also found that variants in CIDEB, which encodes a structural protein found in hepatic lipid droplets, had a protective effect. The burden of rare predicted loss-of-function variants plus missense variants in CIDEB (combined carrier frequency, 0.7%) was associated with decreased alanine aminotransferase levels (beta per allele, -1.24 U per liter; 95% confidence interval [CI], -1.66 to -0.83; P = 4.8×10-9) and with 33% lower odds of liver disease of any cause (odds ratio per allele, 0.67; 95% CI, 0.57 to 0.79; P = 9.9×10-7). Rare coding variants in CIDEB were associated with a decreased risk of liver disease across different underlying causes and different degrees of severity, including cirrhosis of any cause (odds ratio per allele, 0.50; 95% CI, 0.36 to 0.70). Among 3599 patients who had undergone bariatric surgery, rare coding variants in CIDEB were associated with a decreased nonalcoholic fatty liver disease activity score (beta per allele in score units, -0.98; 95% CI, -1.54 to -0.41 [scores range from 0 to 8, with higher scores indicating more severe disease]). In human hepatoma cell lines challenged with oleate, CIDEB small interfering RNA knockdown prevented the buildup of large lipid droplets. CONCLUSIONS: Rare germline mutations in CIDEB conferred substantial protection from liver disease. (Funded by Regeneron Pharmaceuticals.).

UGT1A1 genetic variants are associated with increases in bilirubin levels in rheumatoid arthritis patients treated with sarilumab.

Sarilumab is a human monoclonal antibody against interleukin (IL)-6Rα that has been approved for the treatment of adult patients with moderately to severely active rheumatoid arthritis (RA) and an inadequate response or intolerance to one or more disease-modifying antirheumatic drugs (DMARDs). Mild liver function test abnormalities have been observed in patients treated with sarilumab. We describe a genome-wide association study of bilirubin elevations in RA patients treated with sarilumab. Array genotyping and exome sequencing were performed on DNA samples from 1075 patients. Variants in the UGT1A1 gene were strongly associated with maximum bilirubin elevations in sarilumab-treated patients (rs4148325; p = 2.88 × 10-41) but were not associated with aminotransferase elevations. No other independent loci showed evidence of association with bilirubin elevations after sarilumab treatment. These findings suggest that most bilirubin increases during sarilumab treatment are related to genetic variation in UGT1A1 rather than underlying liver injury.

Clinical Implications of the Amyloidogenic V122I Transthyretin Variant in the General Population.

BACKGROUND: The V122I variant in transthyretin (TTR) is the most common amyloidogenic mutation worldwide. The aim of this study is to describe the cardiac phenotype and risk for adverse cardiovascular outcomes of young V122I TTR carriers in the general population. METHODS AND RESULTS: TTR genotypes were extracted from whole-exome sequence data in participants of the Dallas Heart Study. Participants with African ancestry, available V122I TTR genotypes (N = 1818) and either cardiac magnetic resonance imaging (n = 1364) or long-term follow-up (n = 1532) were included. The prevalence of V122I TTR carriers (45 ± 10 years) was 3.2% (n/N = 59/1818). The V122I TTR carriers had higher baseline left ventricular wall thickness (8.52 ± 1.82 vs 8.21 ± 1.62 mm, adjusted P = .038) than noncarriers, but no differences in other cardiac magnetic resonance imaging measures (P > .05 for all). Although carrier status was not associated with amino terminal pro-B-type natriuretic peptide (NT-proBNP) at baseline (P = .79), V122I TTR carriers had a greater increase in NT-proBNP on follow-up than noncarriers (median 28.5 pg/mL, interquartile range 11.4-104.1 pg/mL vs median 15.9 pg/mL, interquartile range 0.0-43.0 pg/mL, adjusted P = .018). V122I TTR carriers were at a higher adjusted risk of heart failure (hazard ratio 3.82, 95% confidence interval 1.80-8.13, P < .001), cardiovascular death (hazard ratio 2.65, 95% confidence interval 1.14-6.15, P = .023), and all-cause mortality (hazard ratio 1.95, 95% confidence interval 1.08-3.51, P = .026) in comparison with noncarriers. CONCLUSIONS: V122I TTR carrier status was associated with a greater increase in NT-proBNP, slightly greater left ventricular wall thickness, and a higher risk for heart failure, cardiovascular death, and all-cause mortality. These findings suggest the need to develop amyloidosis screening strategies for V122I TTR carriers.

Medical manifestations and health care utilization among adult MyCode participants with neurodevelopmental psychiatric copy number variants.

PURPOSE: Recurrent pathogenic copy number variants (pCNVs) have large-effect impacts on brain function and represent important etiologies of neurodevelopmental psychiatric disorders (NPDs), including autism and schizophrenia. Patterns of health care utilization in adults with pCNVs have gone largely unstudied and are likely to differ in significant ways from those of children. METHODS: We compared the prevalence of NPDs and electronic health record-based medical conditions in 928 adults with 26 pCNVs to a demographically-matched cohort of pCNV-negative controls from >135,000 patient-participants in Geisinger's MyCode Community Health Initiative. We also evaluated 3 quantitative health care utilization measures (outpatient, inpatient, and emergency department visits) in both groups. RESULTS: Adults with pCNVs (24.9%) were more likely than controls (16.0%) to have a documented NPD. They had significantly higher rates of several chronic diseases, including diabetes (29.3% in participants with pCNVs vs 20.4% in participants without pCNVs) and dementia (2.2% in participants with pCNVs vs 1.0% participants without pCNVs), and twice as many annual emergency department visits. CONCLUSION: These findings highlight the potential for genetic information-specifically, pCNVs-to inform the study of health care outcomes and utilization in adults. If, as our findings suggest, adults with pCNVs have poorer health and require disproportionate health care resources, early genetic diagnosis paired with patient-centered interventions may help to anticipate problems, improve outcomes, and reduce the associated economic burden.

Exome sequencing and analysis of 454,787 UK Biobank participants.

A major goal in human genetics is to use natural variation to understand the phenotypic consequences of altering each protein-coding gene in the genome. Here we used exome sequencing1 to explore protein-altering variants and their consequences in 454,787 participants in the UK Biobank study2. We identified 12 million coding variants, including around 1 million loss-of-function and around 1.8 million deleterious missense variants. When these were tested for association with 3,994 health-related traits, we found 564 genes with trait associations at P ≤ 2.18 × 10-11. Rare variant associations were enriched in loci from genome-wide association studies (GWAS), but most (91%) were independent of common variant signals. We discovered several risk-increasing associations with traits related to liver disease, eye disease and cancer, among others, as well as risk-lowering associations for hypertension (SLC9A3R2), diabetes (MAP3K15, FAM234A) and asthma (SLC27A3). Six genes were associated with brain imaging phenotypes, including two involved in neural development (GBE1, PLD1). Of the signals available and powered for replication in an independent cohort, 81% were confirmed; furthermore, association signals were generally consistent across individuals of European, Asian and African ancestry. We illustrate the ability of exome sequencing to identify gene-trait associations, elucidate gene function and pinpoint effector genes that underlie GWAS signals at scale.

Sequencing of 640,000 exomes identifies GPR75 variants associated with protection from obesity.

Large-scale human exome sequencing can identify rare protein-coding variants with a large impact on complex traits such as body adiposity. We sequenced the exomes of 645,626 individuals from the United Kingdom, the United States, and Mexico and estimated associations of rare coding variants with body mass index (BMI). We identified 16 genes with an exome-wide significant association with BMI, including those encoding five brain-expressed G protein-coupled receptors (CALCR, MC4R, GIPR, GPR151, and GPR75). Protein-truncating variants in GPR75 were observed in ~4/10,000 sequenced individuals and were associated with 1.8 kilograms per square meter lower BMI and 54% lower odds of obesity in the heterozygous state. Knock out of Gpr75 in mice resulted in resistance to weight gain and improved glycemic control in a high-fat diet model. Inhibition of GPR75 may provide a therapeutic strategy for obesity.

Acute distal biceps tendon repair using cortical button fixation results in excellent short- and long-term outcomes : a single-centre experience of 102 patients.

AIMS: Acute distal biceps tendon repair reduces fatigue-related pain and minimizes loss of supination of the forearm and strength of flexion of the elbow. We report the short- and long-term outcome following repair using fixation with a cortical button techqniue. METHODS: Between October 2010 and July 2018, 102 patients with a mean age of 43 years (19 to 67), including 101 males, underwent distal biceps tendon repair less than six weeks after the injury, using cortical button fixation. The primary short-term outcome measure was the rate of complications. The primary long-term outcome measure was the abbreviated Disabilities of the Arm, Shoulder and Hand (QuickDASH) score. Secondary outcomes included the Oxford Elbow Score (OES), EuroQol five-dimension three-level score (EQ-5D-3L), satisfaction, and return to function. RESULTS: Eight patients (7.8%) had a major complication and 34 (33.3%) had a minor complication. Major complications included re-rupture (n = 3; 2.9%), unrecovered nerve injury (n = 4; 3.9%), and surgery for heterotopic ossification (n = 1; 1.0%). Three patients (2.9%) overall required further surgery for a complication. Minor complications included neurapraxia (n = 27; 26.5%) and superficial infection (n = 7; 6.9%). A total of 33 nerve injuries occurred in 31 patients (30.4%). At a mean follow-up of five years (1 to 9.8) outcomes were available for 86 patients (84.3%). The median QuickDASH, OES, EQ-5D-3L, and satisfaction scores were 1.2 (IQR 0 to 5.1), 48 (IQR 46 to 48), 0.80 (IQR 0.72 to 1.0), and 100/100 (IQR 90 to 100), respectively. Most patients were able to return to work (81/83, 97.6%) and sport (51/62,82.3%). Unrecovered nerve injury was associated with an inferior outcome according to the QuickDASH (p = 0.005), OES (p = 0.004), EQ-5D-3L (p = 0.010), and satisfaction (p = 0.024). Multiple linear regression analysis identified an unrecovered nerve injury to be strongly associated with an inferior outcome according to the QuickDASH score (p < 0.001), along with infection (p < 0.001), although re-rupture (p = 0.440) and further surgery (p = 0.652) were not. CONCLUSION: Acute distal biceps tendon repair using cortical button fixation was found to result in excellent patient-reported outcomes and health-related quality of life. Although rare, unrecovered nerve injury adversely affects outcome. Cite this article: Bone Joint J 2021;103-B(7):1284-1291.

Advancing human genetics research and drug discovery through exome sequencing of the UK Biobank.

The UK Biobank Exome Sequencing Consortium (UKB-ESC) is a private-public partnership between the UK Biobank (UKB) and eight biopharmaceutical companies that will complete the sequencing of exomes for all ~500,000 UKB participants. Here, we describe the early results from ~200,000 UKB participants and the features of this project that enabled its success. The biopharmaceutical industry has increasingly used human genetics to improve success in drug discovery. Recognizing the need for large-scale human genetics data, as well as the unique value of the data access and contribution terms of the UKB, the UKB-ESC was formed. As a result, exome data from 200,643 UKB enrollees are now available. These data include ~10 million exonic variants-a rich resource of rare coding variation that is particularly valuable for drug discovery. The UKB-ESC precompetitive collaboration has further strengthened academic and industry ties and has provided teams with an opportunity to interact with and learn from the wider research community.

Genome-wide association analysis of serum alanine and aspartate aminotransferase, and the modifying effects of BMI in 388k European individuals.

Serum alanine aminotransferase (ALT) and aspartate aminotransferase (AST) are biomarkers for liver health. Here we report the largest genome-wide association analysis to date of serum ALT and AST levels in over 388k people of European ancestry from UK biobank and DiscovEHR. Eleven million imputed markers with a minor allele frequency (MAF) ≥ 0.5% were analyzed. Overall, 300 ALT and 336 AST independent genome-wide significant associations were identified. Among them, 81 ALT and 61 AST associations are reported for the first time. Genome-wide interaction study identified 9 ALT and 12 AST independent associations significantly modified by body mass index (BMI), including several previously reported potential liver disease therapeutic targets, for example, PNPLA3, HSD17B13, and MARC1. While further work is necessary to understand the effect of ALT and AST-associated variants on liver disease, the weighted burden of significant BMI-modified signals is significantly associated with liver disease outcomes. In summary, this study identifies genetic associations which offer an important step forward in understanding the genetic architecture of serum ALT and AST levels. Significant interactions between BMI and genetic loci not only highlight the important role of adiposity in liver damage but also shed light on the genetic etiology of liver disease in obese individuals.

Computationally efficient whole-genome regression for quantitative and binary traits.

Genome-wide association analysis of cohorts with thousands of phenotypes is computationally expensive, particularly when accounting for sample relatedness or population structure. Here we present a novel machine-learning method called REGENIE for fitting a whole-genome regression model for quantitative and binary phenotypes that is substantially faster than alternatives in multi-trait analyses while maintaining statistical efficiency. The method naturally accommodates parallel analysis of multiple phenotypes and requires only local segments of the genotype matrix to be loaded in memory, in contrast to existing alternatives, which must load genome-wide matrices into memory. This results in substantial savings in compute time and memory usage. We introduce a fast, approximate Firth logistic regression test for unbalanced case-control phenotypes. The method is ideally suited to take advantage of distributed computing frameworks. We demonstrate the accuracy and computational benefits of this approach using the UK Biobank dataset with up to 407,746 individuals.

Mutation spectrum of NOD2 reveals recessive inheritance as a main driver of Early Onset Crohn's Disease.

Inflammatory bowel disease (IBD), clinically defined as Crohn's disease (CD), ulcerative colitis (UC), or IBD-unclassified, results in chronic inflammation of the gastrointestinal tract in genetically susceptible hosts. Pediatric onset IBD represents ≥ 25% of all IBD diagnoses and often presents with intestinal stricturing, perianal disease, and failed response to conventional treatments. NOD2 was the first and is the most replicated locus associated with adult IBD, to date. However, its role in pediatric onset IBD is not well understood. We performed whole-exome sequencing on a cohort of 1,183 patients with pediatric onset IBD (ages 0-18.5 years). We identified 92 probands with biallelic rare and low frequency NOD2 variants accounting for approximately 8% of our cohort, suggesting a Mendelian inheritance pattern of disease. Additionally, we investigated the contribution of recessive inheritance of NOD2 alleles in adult IBD patients from a large clinical population cohort. We found that recessive inheritance of NOD2 variants explains ~ 7% of cases in this adult IBD cohort, including ~ 10% of CD cases, confirming the observations from our pediatric IBD cohort. Exploration of EHR data showed that several of these adult IBD patients obtained their initial IBD diagnosis before 18 years of age, consistent with early onset disease. While it has been previously reported that carriers of more than one NOD2 risk alleles have increased susceptibility to Crohn's Disease (CD), our data formally demonstrate that recessive inheritance of NOD2 alleles is a mechanistic driver of early onset IBD, specifically CD, likely due to loss of NOD2 protein function. Collectively, our findings show that recessive inheritance of rare and low frequency deleterious NOD2 variants account for 7-10% of CD cases and implicate NOD2 as a Mendelian disease gene for early onset Crohn's Disease.

Heterozygosity for a Pathogenic Variant in SLC12A3 That Causes Autosomal Recessive Gitelman Syndrome Is Associated with Lower Serum Potassium.

BACKGROUND: Potassium levels regulate multiple physiologic processes. The heritability of serum potassium level is moderate, with published estimates varying from 17% to 60%, suggesting genetic influences. However, the genetic determinants of potassium levels are not generally known. METHODS: A whole-exome sequencing association study of serum potassium levels in 5812 subjects of the Old Order Amish was performed. A dietary salt intervention in 533 Amish subjects estimated interaction between p.R642G and sodium intake. RESULTS: A cluster of variants, spanning approximately 537 kb on chromosome 16q13, was significantly associated with serum potassium levels. Among the associated variants, a known pathogenic variant of autosomal recessive Gitelman syndrome (p.R642G SLC12A3) was most likely causal; there were no homozygotes in our sample. Heterozygosity for p.R642G was also associated with lower chloride levels, but not with sodium levels. Notably, p.R642G showed a novel association with lower serum BUN levels. Heterozygotes for p.R642G had a two-fold higher rate of self-reported bone fractures and had higher resting heart rates on a low-salt diet compared with noncarriers. CONCLUSIONS: This study provides evidence that heterozygosity for a pathogenic variant in SLC12A3 causing Gitelman syndrome, a canonically recessive disorder, contributes to serum potassium concentration. The findings provide insights into SLC12A3 biology and the effects of heterozygosity on electrolyte homeostasis and related subclinical phenotypes that may have implications for personalized medicine and nutrition.

Peat-forest burning smoke in Maritime Continent: Impacts on receptor PM2.5 and implications at emission sources.

This study characterizes the impacts of transported peat-forest (PF) burning smoke on an urban environment and evaluates associated source burning conditions based on carbon properties of PM2.5 at the receptor site. We developed and validated a three-step classification that enables systematic and more rapid identification of PF smoke impacts on a tropical urban environment with diverse emissions and complex atmospheric processes. This approach was used to characterize over 300 daily PM2.5 data collected during 2011-2013, 2015 and 2019 in Singapore. A levoglucosan concentration of ≥0.1 µg/m3 criterion indicates dominant impacts of transported PF smoke on urban fine aerosols. This approach can be used in other ambient environments for practical and location-dependent applications. Organic carbon (OC) concentrations (as OC indicator) can be an alternate to levoglucosan for assessing smoke impacts on urban environments. Applying the OC concentration indicator identifies smoke impacts on ∼80% of daily samples in 2019 and shows an accuracy of 51-86% for hourly evaluation. Following the systematic identification of urban PM2.5 predominantly affected by PF smoke in 2011-2013, 2015 and 2019, we assessed the concentration ratio of char-EC/soot-EC as an indicator of smoldering- or flaming-dominated burning emissions. When under the influence of transported PF smoke, the mean concentration ratio of char-EC to soot-EC in urban PM2.5 decreased by >70% from 8.2 in 2011 to 2.3 in 2015 but increased to 3.8 in 2019 (p < 0.05). The reversed trend with a 65% increase from 2015 to 2019 shows stronger smoldering relative to flaming, indicating a higher level of soil moisture at smoke origins, possibly associated with rewetting and revegetating peatlands since 2016.